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Abstract 

Assuming two-neutrino mixing takes place in vacuum, we study in detail the 
conditions under which the v<i — ► v e and — ► v e (y e — ► ^( T )) transitions in 
the Earth are strongly enhanced by the neutrino oscillation length resonance 
when the neutrinos cross the Earth core. We show, in particular, that the 
neutrino oscillation length resonance is operative also in the — > D s (or 
ffj. v s) transitions at small mixing angles. The properties of the v?, — * 
v e and Up — > v e (y e — > v^m) transition probabilities in the corresponding 
resonance regions are examined. Some implications of our results for the 
transitions of solar and atmospheric neutrinos traversing the Earth core, 
relevant for the interpretation of the results of the solar and atmospheric 
neutrino experiments, are also discussed. 



*Also at: Institute of Nuclear Research and Nuclear Energy, Bulgarian Academy of Sciences, 1784 
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1. Introduction 

The Vi — > v e and — + v e [y e — > ^ M ( T )) transitions/oscillations of solar and atmospheric 
neutrinos in the Earth, caused by neutrino mixing (with nonzero mass neutrinos) in vacuum 
Q can be strongly amplified by a new type of resonance which differs from the MSW one and 
takes place when the neutrinos traverse the Earth core on the way to the detector [|]]. At 
small mixing angles (sin 2 29 ^ 0.05), the maxima due to this resonance in the corresponding 
transition probabilities, P{y 2 — > v e ) = P e 2 and Piy^ — > u e ) [P{y e — * ^(t)))> are absolute 
maxima and dominate in P e2 and Piv^ —>■ u e ): the values of the probabilities at these 
maxima in the simplest case of two-neutrino mixing are considerably larger - by a factor of 
~ (2.5 - 4.0) (~ (3.0 - 7.0)), than the values of P e2 and Pfa -> u e ) = P{v e -> z/ m(t) ) at 
the local maxima associated with the MSW effect taking place in the Earth core (mantle). 
The enhancement is less dramatic at large mixing angles. Even at small mixing angles the 
resonance is relatively wide in the neutrino energy (or resonance density) - it is somewhat 
wider than the MSW resonance. It also exhibits strong energy dependence. 

The conditions for existence of the indicated resonance in the probabilities P e2 and 
P(z/ M — > u e ) (P(v e — > Vn( T ))) include specific constraints on the neutrino oscillation length 
in the Earth mantle and in the Earth core JTJ. When satisfied, these conditions ensure that 
the relevant oscillating factors in the v 2 — > v e and — > v e [y e —>■ u^ T )) transition proba- 
bilities are maximal and that this produces a resonance maximum in P e2 and Piy^ — > v e ) 
{P{v e ~ > u h(t)))- Accordingly, the term "neutrino oscillation length resonance" or simply 
"oscillation length resonance" was used in |TJ to denote the resonance of interest In con- 
trast, the MSW effect is, as is well-known, a resonance in the neutrino mixing. There exists 
a beautiful analogy between the neutrino oscillation length resonance and the electron spin- 
flip resonance taking place in a specific configuration of two magnetic fields (see for a 
somewhat more detailed discussion) 0. The neutrino oscillation length resonance is also in 
many respects similar to the resonance transforming a circularly polarized LH photon into 
a RH photon when the photon traverses three layers of optically active medium such that 
the optical activity and the length of the path of the photon in the first and the third layers 
are identical but differ from those in the second layer. 

The implications of the oscillation length resonance enhancement of the probability P e2 
for the tests of the MSW v e — > v^u) and v e v s transition solutions of the solar neutrino 
problem are discussed in ref. and in much greater detail in refs. (see also 0]). It is 
quite remarkable that for values of Am 2 = (4.0 — 8.0) x 10 -6 eV 2 from the small mixing angle 
(SMA) MSW solution region and the geographical latitudes at which the Super-Kamiokande, 
SNO and ICARUS detectors are located, the enhancement takes place in the v e — > v^{ T ) case 
for values of the 8 B neutrino energy lying in the interval ~ (5 — 12) MeV to which these 



x The u 2 — * v e transition probability accounts, as is well-known, for the Earth effect in the solar 
neutrino survival probability in the case of the MSW two- neutrino v e — * zv T ) and v e — * v s transition 
solutions of the solar neutrino problem, v s being a sterile neutrino. 

2 The oscillation length resonance can produce deep absolute minima in the probability P e2 as well 
(see further). 

3 This analogy was brought to the attention of the author of [Q] by L. Wolfenstein. 
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detectors are sensitive. The resonance maximum in P e 2 at sin 2 29 = 0.01 for the trajectory 
with a Nadir angle h = 23°, for instance, is located at a resonance density pj^an ~ 8-0 g/cm 3 , 
which corresponds to E = 5.3 (10.5) MeV if Am 2 = 4.0 (8.0) x 10" 6 eV 2 . Accordingly, 
at small mixing angles this enhancement is predicted to produce [@] a much bigger - by a 
factor of ~ 6, day- night (D-N) asymmetry in the Super-Kamiokande sample of solar neutrino 
events, whose night fraction is due to the core-crossing solar neutrinos, in comparison with 
the asymmetry determined by using the whole night Super-Kamiokande event sample On 
the basis of these results it was concluded in that it may be possible to test a rather 
large part of the MSW v e — * u^m small mixing angle (SMA) solution region in the Am 2 — 
sin 2 26 l plane by performing selective, i.e., core D-N asymmetry measurements. The Super- 
Kamiokande collaboration has already successfully applied this approach to the analysis of 
their solar neutrino data in terms of the MSW effect |||7]]: the limit the collaboration has 
obtained on the D-N asymmetry utilizing only the core event sample permitted to exclude 
a part of the SMA solution region allowed by the mean event rate data from all solar 
neutrino experiments (Homestake, GALLEX, SAGE, Kamiokande and Super-Kamiokande). 
The SMA solution values of sin 2 29 and Am 2 thus excluded belong to the intervals sin 2 26 = 
(0.007 — 0.01) and Am 2 = (5.0 x 10~ 6 — 10 -5 ) eV 2 . In contrast, the current upper limit on 
the D-N asymmetry, derived by the Super-Kamiokande collaboration using the whole night 
data sample of solar neutrino events, does not permit to probe the SMA solution region: 
the predicted asymmetry is too small @. The oscillation length resonance leads also to a 
relatively large D-N asymmetry in the recoil - e~ spectrum which is being measured in the 
Super-Kamiokande experiment @. 

As was shown in [|1], the conditions for the neutrino oscillation length resonance are not 
even approximately fulfilled if solar neutrinos take part in v e — > v s transitions, and hence the 
effect of strong enhancement of the probability P e 2 does not take place in this case. However, 
when the MSW resonance occurs in the Earth core, the purely MSW enhancement of P e 2 
at small mixing angles is rather strongly "assisted" by some of the additional interference 
terms present in P e 2 (see ref. [0 for details) and at the corresponding maximum P e 2 has a 
noticeably larger value (by a factor of ~ (2 — 4)) than if it were determined by the MSW 
effect only. Still, at sin 2 29 ^ 0.02, the value of P e 2 at its absolute maximum in the v e — > v s 
case is substantially smaller (by a factor of ~ (2.5 — 4.0) [|],[§) than the maximum value of 
P e 2 due to the oscillation length resonance if v e — > v^r) transitions took place. In addition, 
due basically to the well-known difference between the v^m and u s effective potentials in 
matter and the fact that the Earth matter is approximately isotopically symmetric, for fixed 
Am 2 and sin 2 29 the absolute maximum of P e 2 occurs in the case of v e — > v s transitions at 
values of the neutrino energy E exceeding approximately by a factor 2 the values at which 
the oscillation length resonance maximum takes place if the solar v e undergo transitions 



4 The possibility of enhancement of the Earth effect in the Super-Kamiokande sample of events 
generated by the core-crossing solar neutrinos was discussed in ref. ||, where on the basis of 
results obtained for few pairs of values of Am 2 and sin 2 29 it was also suggested that the effect 
might be measurable even in the case of small neutrino mixing angle. However, no estimates of the 
magnitude of the core enhancement were given in this article and the enhancement was interpreted 
to be due to the MSW effect taking place in the Earth core. 
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into Vfi{ T ) @- As a consequence, the oscillation length resonance does not produce such a 
dramatic enhancement of the D-N asymmetry in the core sample of the Super-Kamiokande 
solar neutrino events as in the case of v e — > v^r) transitions. Thus, the possibility to probe 
the MSW v e — > v s (SMA) solution of the solar neutrino problem via the measurement of 
the D-N effect related observables in the Super-Kamiokande experiment appears unlikely at 
present ||. 

The oscillation length resonance enhancement can be present also in the —>■ v e and 
v e — > v^t) transitions of atmospheric neutrinos crossing the Earth core |l| . It can take place 
practically for all neutrino trajectories through the core, in particular, for the trajectories 
with h = (0° — 23°). In the case of two-neutrino mixing and at small mixing angles, the 
maximum in Pivp, — > v e ) = P{ye ~* ^(r)) due to the oscillation length resonance is an 
absolute maximum and dominates in these probabilities []. This can be especially relevant for 
the interpretation of the results of the atmospheric neutrino experiments and for the future 
studies of the oscillations/transitions of atmospheric neutrinos crossing the Earth. As is 
well-known, the Super-Kamiokande collaboration has reported recently strong experimental 
evidences for oscillations of the atmospheric (z/^) neutrinos Assuming two-neutrino 
mixing, the data is best described in terms of <-> v T and <-» v T vacuum oscillations with 
parameters Am 2 = (5x 1CT 4 — 6x 1CT 3 ) eV 2 and sin 2 29 = (0.8 — 1.0), although the hypothesis 
of large mixing <-» v s oscillations provides almost equally good fit to the data. At the 
same time, the possibility of two-neutrino v e and <-> u e large mixing oscillations is 

disfavored by the presently existing Super-Kamiokande data |8|||. At Am 2 ^2x 10~ 3 eV 2 
this possibility is ruled out by the results of the CHOOZ experiment [[UJ. 

It is a remarkable coincidence that, as was shown in M, for values of Am 2 ~ (5 x 
10~ 4 — 5 x 10~ 3 ) eV 2 and small mixing angles, sin 2 29 ^ 0.10, the neutrino oscillation length 
resonance in P{y^ — > u e ) (P(v e — ¥ v h{t))) occurs for values of the energy E of the atmospheric 
u e and which contribute either to the sub-GeV or to the multi-GeV samples of e— like 
and fi— like events in the Super-Kamiokande experiment. For sin 2 29 = 0.01, Am 2 = 5 x 
10~ 4 ; 10 -3 ; 5 x 10~ 3 eV 2 , and h — 0° (Earth center crossing), for instance, the absolute 
maximum in Piv^ — > u e ) = P{y e — > ^(t)) due to the oscillation length resonance takes 
place at E = 0.75; 1.50; 7.5 GeV. Thus, for values of Am 2 ~(5x 10~ 4 - 6 x lO" 3 ) eV 2 
of the region of the <-» v T oscillation solution of the atmospheric neutrino problem, the 
oscillation length resonance strongly enhances the — > v e (and v e — > ^W)) transitions 
of the atmospheric neutrinos crossing the Earth core, making the transition probabilities 
large and the transitions perhaps detectable even at small mixing angles. Actually, it was 
suggested in [0,|ll|] that the small excess of e-like events in the region —1 < cos9 z < —0.6, 
9 Z being the Zenith angle, either in the sub-GeV or in the multi-GeV sample of atmospheric 
neutrino events, observed (in both samples) in the Super-Kamiokande experiment [§,0, is 
due to z/^ — > v e small mixing angle, sin 2 2# eAt = (0.01 — 0.10), transitions with Am 2 ~ 
(5 x 10 -4 — 10~ 3 ) eV 2 or respectively Am 2 ~ (2 x 10~ 3 — 6 x 10~ 3 ) eV 2 , strongly enhanced by 
the neutrino oscillation length resonance as neutrinos cross the Earth core on the way to the 



5 Let us note that the oscillation length resonance enhancement of P{v^ — > v e ) is operative at large 
mixing angles as well, at which the probabilities P e 2 and -P(f^ — ► v e ) exhibit different dependence 
on E/ Am 2 (see further). 
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detector. The same resonantly enhanced transitions with Am 2 ~ (2 x 10~ 3 — 6 x 10~ 3 ) eV 2 
(Am 2 ~ (5 x 10 -4 — 10~ 3 ) eV 2 ) should produce [fT|JTlT| at least part of the strong Zenith angle 
dependence, exhibited by the /x— like multi-GeV (sub-GeV) Super-Kamiokande data. 

The — > u e and u e — *> v^m transitions of atmospheric neutrinos under discussion should 
exist if three-flavour-neutrino (or four-neutrino) mixing takes place in vacuum. More specif- 
ically, such transitions are naturally predicted to take place in the case of three-neutrino 
mixing M if large mixing angle v e <-> oscillations with Am 2l ~ 10~ 10 eV 2 or small mixing 
angle MSW v e — > transitions with Am^ ~ (4 — 8) x 10~ 6 eV 2 , provide the solution of the 
solar neutrino problem, and if the atmospheric neutrino anomaly is due to <-> v T large mix- 
ing angle oscillations with Am 31 ~ (5 x 10~ 4 — 6 x 10 -3 ) eV 2 . If Am 31 ^> Am 21 , which seems 
to be a plausible possibility in view of the neutrino oscillation interpretation of the solar and 
atmospheric neutrino data, the relevant three-neutrino — > v e and v e — > v^u) transition 
probabilities reduce to a two-neutrino transition probability (see, e.g., [0) with Am 31 and 
sin 2 28 13 = sin 2 2^ = 4|£7 e3 | 2 (l — |t7, 3 | 2 ) playing the role of the relevant two-neutrino os- 
cillation parameters, where is the element of the lepton mixing matrix U, coupling the 
electron and the third (heaviest) massive neutrino, u 3 , in the charged lepton current. Un- 
der the above assumptions and taking into account the CHOOZ result JL0| one finds that 
\U e3 \ 2 cannot be large, but is also not necessarily exceedingly small: sin 2 26> 13 ^ 0.20. Thus, 
searching for the — ► u e and v e —>■ Vfj T ) transitions of atmospheric neutrinos, amplified by 
the oscillation length resonance, can provide unique information, in particular, about the 



magnitude of the element U e3 of the lepton mixing matrix [|T^]. As is not difficult to show 
using the results of [Q, the fluxes of atmospheric v e ^ having energy E and crossing the 



Earth along a trajectory with Zenith angle 6 Z before reaching the detector, <& Ve (E, 9 Z ) } are 
given by the following expressions in the three- neutrino mixing scheme under discussion [[Uj : 

$ Ve (E,9 z ) = d>° (l + [s 2 23 r(E,9 z ) - 1] P 2u (E,9 z ;Am 2 31 ,9 13 )) , (la) 



e z ) - & (1 + 4 [(s 2 23 r(E, e z ))-' - 1] p 2v (e, e z - a™^, e 



13; 



-2r 2 <? 2 
zc 23*23 



I -Re (e- lK A 2u (u T ^u T )) ), (lb) 



Here $° = 3>° e( . (E, 6 Z ) is the v e {n) flux in the absence of neutrino transitions/oscillations, 

U^ 3 being the //— 1/ 3 element of {7, c 2 3 = 1— sf 3 , P 2v (E, 9 Z \ Am 31 , 6*i 3 ) is the probability of two- 
neutrino transitions in the Earth which coincides in form with P[y e — > u T ), P 2v = P[y e — > 
z/ T ), but describes z/ e — > v T transitions only in a certain limit |T^|, and k and A 2v [y T — > z/ T ) are 
known phase and two-neutrino transition probability amplitude []. The interpretation of the 
Super-Kamiokande atmospheric neutrino data in terms of <-» v T oscillations requires the 
parameter s 2 3 to lie approximately in the interval (0.30 - 0.70), with 0.5 being the statistically 



6 Analytic expressions for P(v e — ► u T ), k and A 2v (y r — ► u T ) will be given later (see eqs. (13), (15) 
and (16). 
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preferred value. For the predicted ratio r(E, 9 Z ) of the atmospheric and v e fluxes for the 
Earth core crossing neutrinos, —1 < cos# z ^ —0.8, one has [14jl5||: r(E,6 z ) = (2.0 — 2.5) for 
the neutrinos giving the contribution to the sub-GeV samples of Super-Kamiokande events, 
and r(E, 9 Z ) = (2.6 — 4.5) for those giving the main contribution to the multi-GeV samples. 
If s 2 23 = 0.5 and r(E,9 z ) = 2.0 we have (sj 3 r(E,9 z ) - 1) ^ and the possible effects of 
the Ufj, —> v e and v e — > v^t) transitions on the v e flux and correspondingly on the sub-GeV 
e— like sample of events, would be strongly suppressed even if these transitions are maximally 
enhanced by the oscillation length resonance. The indicated suppression may actually be 
taking place []. However, one should keep in mind that the factor (s| 3 r(E,6 z ) — 1) can be 
as large as ~ 0.6 for the sub-GeV neutrinos and that the effects of the resonantly enhanced 
— ► v e and v e — > u^ T ) transitions in the sub-GeV e— like event sample can be suppressed by 
the specific way this sample is selected from the data in the Super-Kamiokande experiment 
13[. For the multi-GeV neutrinos we have (s 2 23 r(E,9 z ) - 1) £ 0.3 (0.9) for s 2 23 = 0.5 (0.7). 



The factor (s| 3 r(E,9 z ) — 1) amplifies the effect of the — > v e transitions in the e— like 



sample for E ^ (5 — 6) GeV, for which r(E, 6 Z ) ^ 4 [0,0. This discussion suggests that the 
effects of the neutrino oscillation length resonance in the Super-Kamiokande e— like multi- 
GeV data for Am 31 ~ (2 x 10~ 3 — 6 x 10~ 3 ) eV 2 can be much larger than those in the 
e— like sub-GeV data for Am^ ~(5x 10~ 4 — 10~ 3 ) eV 2 [ |rT| , |T3"| ] ■ As similar analysis of the 
expression (lb) for shows, they are also expected to be larger (within the three-neutrino 
mixing scheme considered) than the oscillation length resonance effects in the multi-GeV 
ji— like event sample |pT|[13| . 

In view of the possible important phenomenological implications of the neutrino oscilla- 
tion length resonance, some of which have been discussed above, we continue in this article 
the studies of the resonance in the v 2 — > v e and — > v e [y e — > v^t)) transitions, began in ref. 
]T|]. Assuming two-neutrino v e — v^( T ) or v e — v s mixing (with nonzero-mass neutrinos) exists 
in vacuum, we study in greater detail the conditions under which the neutrino oscillation 
length resonance takes place in the v 2 — > v e and v e [v e — > ^( T )) transitions of the 

solar and atmospheric neutrinos in the Earth. We give a graphical representation of both 
the necessary conditions for existence of the extrema in the probabilities P e2 and Piy^ —> v e ) 
(P(i/ e — > ^( r )), associated with the resonance, and of the supplementary conditions ensur- 
ing that these extrema correspond to maxima. This permits, e.g., to determine all relevant 
values of the transition parameters, E, Am 2 and sin 2 29, at which the resonance occurs for 
a given neutrino trajectory through the Earth core. Similar analysis is performed for the 
case of z/ M — > v s transitions at small mixing angles. We also investigate the behavior of the 
probabilities P e2 and Piy^ — > u e ) [P{y e v h{t)) i n the region of the neutrino oscillation 
length resonance. We show, in particular, that at small mixing angles the resonance is rel- 
atively wide not only in the variable E/Am 2 , but also in the Nadir angle variable h. The 
implications of our results for the current and future experiments with solar and atmospheric 
neutrinos are also briefly discussed. 



7 Indeed, a more detailed investigation performed within the indicated three-neutrino mixing 
scheme reveals, in particular, that the excess of e-like events in the Super-Kamiokande sub-GeV 
data at —1 < cos6 z < —0.6 seems unlikely to be due to small mixing angle, sin 2 2#i3 ^ 0.20, 
v a v e transitions amplified by the oscillation length resonance. 
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Let us note that the Earth enhancement of the two-neutrino v 2 — > v e and zv e ) — > v e U) 
{y e — > ^(t)) transitions of interest of the solar and atmospheric neutrinos in the case of 
relatively small mixing angles has been discussed in the past in, e.g., refs. PJ5|]T6|-|T9| . Some 



of these articles contain plots of the probabilities P e2 and/or P{y u — > v e ) or P{y e — > ^( r )) on 
which one can clearly recognize now the dominating maximum due to the neutrino oscillation 



length resonance (see, e.g., |]16|-|19|). However, this maximum was invariably interpreted to 



be due to the MSW effect in the Earth core before the appearance of the study |TJ. 

2. The Earth Model, the Two-Layer Density Approximation and 
the Neutrino Oscillation Length Resonance 

We begin with a brief summary of the main properties of the reference Earth model and 
of the results obtained in ref. PJ, which will be utilized in our further analysis. Following 
|I|-|4|,|T7|| we will use the Stacey model from 1977 [20] as a reference Earth model. The density 



distribution in the Stacey model is, as in all Earth models known to us, spherically symmetric 
and there are two major density structures - the core and the mantle, and a certain number 
of substructures (shells or layers). The core has a radius R c = 3485.7 km, the Earth mantle 
depth is approximately R ma n = 2885.3 km, and the Earth radius in the Stacey model is 
i?9 = 6371 km. The mean values of the matter densities of the core and of the mantle read, 
respectively: p c = 11.5 g/cm 3 and p man = 4.5 g/cm 3 . The density distribution in the 1977 



Stacey model practically coincides with that in the more recent PREM model |21] . 

All the interesting features of the solar and atmospheric neutrino transitions in the Earth, 
including those related to the neutrino oscillation length resonance, can be understood quan- 
titatively in the framework of the two-layer model of the Earth density distribution The 
density profile of the Earth in the two-layer model is assumed to consist of two structures 
- the mantle and the core, having different densities, p man and p c , and different electron 
fraction numbers, Y e man and Y" e c , none of which however vary within a given structure. The 
core radius and the depth of the mantle are known with a rather good precision and these 



data are incorporated in the Earth models pOJSIJ. The densities p man and p c in the case 



of interest should be considered as mean effective densities along the neutrino trajectories, 
which can vary somewhat with the change of the trajectory p man = p m an and p c = p c . In 
the Stacey model one has: p man = (4 — 5) g/cm 3 and p c = (11 — 12) g/cm 3 . For the electron 
fraction numbers in the mantle and in the core one can use the standard values [F2T1H2~I3] (see 



also @) Y e man = 0.49 and Y e c = 0.467. Numerical calculations show |3| (see also [0) that, 
e.g., the (time-averaged probability P e2 calculated within the two- layer model of the Earth 
with Pman and p c taken from the Stacey 1977 model |2(J reproduces with a remarkably high 
precision the probability calculated by solving numerically the relevant system of evolution 
equations with the much more sophisticated Earth density profile of the Stacey model 



The probability of the v 2 — > v e transition in the Earth, P e2 , which accounts for the Earth 
effect in the solar neutrino survival probability, has the following form in the two-layer model 



8 Let us note that because of the spherical symmetry of the Earth, a given neutrino trajectory 
through the Earth is completely specified by its Nadir angle (see, e.g., §). 
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of the Earth density distribution ffl: 



1 



P e2 = sin 2 9 + - [1 - cos A£"X"] [sin 2 (20, 



m 



sin 2 # 



+ - [1 — cos A£"X"] [1 - cos A£'X'] [sin 2 (2C - A9' m + 9)- sin 2 (20. 



2/na" 

rn 



1 



[1 - cos AE"X"] [1 - cos2A£X'] sin 2 (2^ - 9) - sin 2 cos 2 (2C - 26' m ) 



+ - [1 + cos AE"X"] [1 - cos2A£'X'] [sin 2 (2^ 



+ - sin AE"X" sin 2 AE'X' 
2 



sm 2 (2^ 



sin 2 fl 



1) - sin 2 
cos(2C - 29' m 



+ - [cos(A£'X' - AE"X") - cos(AE'X' + AE"X")] sm(49' m -29) sm(29'^-29' n 



Here 



AE' (AE h 



Am 2 



2E \ 



1 



Pman (c) 



cos 2 29 + sin 2 29 , 



(3) 



Pman (c) , 

is the difference between the energies of the two energy- eigenstate neutrinos in the mantle 
(core), 9' m and 9'^ are the mixing angles in matter in the mantle and in the core, respectively, 



sin 2 29L (sin 2 29. 



2 c-.nl I 



sin 2 29 



P manic) \ 2 2 nn , ■ 2 nn 

cos 29 + sin 29 



(4) 



'1 - 



Pman(c) 

X' is half of the distance the neutrino travels in the mantle and X" is the length of the path 
of the neutrino in the core, and p r ^ an and p r c es are the resonance densities in the mantle and 
in the core. The latter can be obtained from the expressions 



P 



Am 2 cos 29 
2EV2G F f(Y e ) 



m N , 



f(Y e ) = Y e , u e — > u^ T ) {.v u -> u e ) transitions, 



/(F e ) = -(3F e -l), u e 



u s \v R 



u P ) transitions, 



(5) 

(6a) 
(66) 



m^r being the nucleon mass, by using the specific values of Y e in the mantle and in the core. 
We have p™ s an ^ p r c es (but p^ s an ~ p r c es ) because Y e c = 0.467 and Y e man = 0.49. Obviously, 
one has: f(Y e c )p^ = f(Yr n )Pman- 

We will present some of our results (obtained for fixed sin 2 29) utilizing the density 
parameter 

Am 2 COS26 1 



Pr 



2EV2G f C, 



-m N 



(7a) 



F^a's) 



where the constant C a i s \ = 0.50 (0.25) for the u e — > v^m {v e - > v s ) transitions. The 
parameter p r would coincide with the resonance density if Y e were equal to 1/2 both in the 
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mantle and in the core, so p r = p^lni f° r given sin 2 29 it is equivalent to E / Am 2 , but gives 
an idea about the densities at which one has an enhancement of the transition probabilities 
of interest when the neutrinos cross the Earth and we have used it for the latter purpose 
in |H and in @. Clearly, C aPr = Y e c p r c es = Y e man p™ s an , or C s p r = 0.5(3y e c - l)p r c es = 
0.5(3F e man — l)Pman depending on whether v e — v^( T ) or v e — v s mixing takes place. The 
relation between p r and Am 2 /E, expressed in units of g/cm 3 and eV 2 /MeV, respectively, 
has the form: 

For a neutrino trajectory which is specified by a given Nadir angle h we have: 

X' = R e cosh -^R 2 C -R% sin 2 h, X" = 2 ^ R 2 C - R% sin 2 h. (8) 

The conditions for the presence of a neutrino oscillation length resonance maximum in 
Pe2 read 0: 

AE'X' = n(2k + 1), AE"X" = n(2k' + 1), k, k' = 0, 1, 2, ... , (9) 

and 

sin 2 (2^ - A9' m + 9)- sin 2 9 > 0, (10a) 
sin 2 (2C - + 9) sin(2C - 29'J sm(29' rn - 9) cos(2^ - A9' m + 9) 

+ - sin 2 9 sin(4C - 80^ + 29) sin(40£ - < 0. (106) 

When equalities (9) hold, conditions (10a) and (10b) ensure that P e2 has a maximum 0. At 
the resonance maximum the probability P e2 takes the value [TJ: 

= sin 2 (2C-4C + ^)- (11) 
If instead of inequality (10a) we have 

sin 2 (2C - A9' m + 9)- sin 2 9 < 0, (12) 

and inequality (10b) is fulfilled (together with equalities (9)), the neutrino oscillation reso- 
nance extrema in P e2 correspond to minima. 

As was already indicated in the Introduction, the oscillation length resonance is present 
in the — > v e (and v e — > i^m) transitions of atmospheric neutrinos as well JT|. The 
probability of the — > i> e (z/ e — > i^m) transitions, P(z^ M — > ^ e ) = P (f e "~ > v v-W))-> can De 
obtained from eq. (1) by formally setting 9 = while keeping 9' m ^ and 6^ 7^ 0: 

P(^ -> ^e) = P(^e -> ^(r)) = £ [1 ~ COS AP"X"] sin 2 2^ 



9 Note that the factor 1/4 in the second term in the left-hand side of the inequality (10b) is missing 
in the same inequality given in [U (see eq. (12b) in |fj]). 

10 All subsequent results for the — ► v e and v e — ► ^( r ) transitions will be formulated for the 
probability Piy^ — > v e ). 
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+ - [1 - cos AE"X"] [1 - cos AE'X'} [sin 2 (2C - AB' m ) - sin 2 29" m 



[1 - cos AE"X"} [1 - cos 2 AE'X'] sin 2 28' m cos 2 (20^ - 26' m 



+ | [1 + cos A£"X"] [1 - cos 2 AE'X'] sin 2 29' m 
+ ^sinA£"X" sin 2 AE'X' sin 2 20^ cos(20^ - 29'J 

+ - [cos{AE'X' - AE"X") - cos(A£'X' + AE"X")} sin A6' m sin(20^ - 26'J. (13) 

Correspondingly, the maximum conditions (9) are the same as for P e2 , while the analogs 
of the conditions (10a) and (10b) for this probability change: condition (10a) becomes 
sin 2 (20^ — 4:9' m ) > and is always fulfilled, while condition (10b) transforms into [0] 



cos(20 



40 < 0. 



(14) 



Obviously, in contrast to P e 2, the probability Piy^ — > z/ e ) cannot have minima associated 
with conditions (9). 

At small mixing angles, sin 2 29 < 0.10, the maxima in the probabilities P e 2 and Piy^ — > 
u e ) associated with the conditions (9) and (10a), (10b), and (9) and (14), are absolute 
maxima. For given sin 2 29 and trajectory crossing the Earth core, the neutrino oscillation 
length resonance is relatively wide in the variable E/Am 2 (or p r , see Figs. 1 and 2 in [l|]). 

The analytic expression and the resonance conditions (13) and (14) derived for the 
probability P[y a — > u e ) are valid, with minor modifications, for the probabilities of the 
i) v e — > v s {p e —>■ u s ) and ii) — > u s [p n — > u s ) transitions, P{y e — * v s ) {P(v e ~> v s )) and 
P{ v ii ~^ v s) {PiPn ~~ * ^s)) : one has to use in eqs. (3), (4), (9) and (14) and in the formula for 
P{ v u ~ * v e)i the expressions for the resonance densities p: 



res 

man,c 



with i) f(Y e ) given by eq. (6) 
(f(Y e ) = 0.5(1 - 3F e )) and ii) f(Y e ) = 0.5(Y e - 1) (f(Y e ) = 0.5(1 - Y e )), respectively. Since 
Y e for the Earth is close to 0.5, only the probabilities P[y e — > v s ) and Piv^ — > u s ) can, in 
principle, be amplified for cos 29 > by the oscillation length resonance if Am 2 > 0, which 
is assumed throughout this study. Note, however, that in the case of two-neutrino mixing, 
the probability P{y e — > u s ) (P{p n — > u s )) for Am 2 > coincides, as can be easily shown, 
with the probability P{p e — ► v s ) (P(f M — > u s )) for Am 2 < for any given values of E, |Am 2 | 
and 9, and that conditions (9) and (14) are identical in the two cases. 

The phase k and the probability amplitude A2 V (v T — > v T ), which appear in eq. (lb) for 
the flux of atmospheric in the case of three-flavour neutrino mixing and strong hierarchy 
between the neutrino mass squared differences and therefore can play important role in the 
interpretation of the, e.g., Super-Kamiokande atmospheric neutrino data, have the following 
form in the two-layer model of the Earth density distribution flI,IT|,r2 : 



« = \\~^ X+V2G F — (X"Y:pc+2X'Y: 



man — \ 
Pman ) 



-2AP'X'-AP"X"]-=^1 
J 2E 



i2AE'X' _ j 



1 + (e-* 



AE"X'' 



1 cos 2 (^-fi£ 



X COS 2<fii2, 

(15) 

2 ni 



cos' 9 
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+ (e- lAB "*"-l)cos 2 C + \ (e- iAE " X "-l) ( e -^'-l)sin(2C-2C)sin2C 

(16) 

where X = X" + 2X' and[]cos2^ 12 = (\U el \ 2 - \U e2 \ 2 )/(l - \U e3 \ 2 ). 

Let us note finally that the neutrino oscillation length resonance arises when the condi- 
tions specified above are fulfilled, due basically to the interference between the probability 
amplitudes describing the neutrino transitions in the Earth core and in the Earth mantle. 

3. The Oscillation Length Resonance Conditions and the 

Probabilities P e2 , P(v a -> u e ) (P(u e -> v^ T ))) and P[p n -> u s ) 

3.1. The Resonance Conditions 

The quantities AE'X' and AE"X" represent the phase differences the states of the 
energy-eigenstate neutrinos acquire after the neutrinos have crossed the mantle and the 
core, respectively. Thus, the neutrino osacillation length resonance can take place only when 
these phases are correlated, being odd multiples of tt, eq. (9). Since 

A# (AE") = — , (17) 

man(c) 

where L™ an i c } is the neutrino oscillation length in matter in the Earth mantle (core), and X' 
and X" are fixed for a given neutrino trajectory, conditions (9) are constraints on the the 
oscillation lengths L™ an and L™. When these constraints are satisfied, the oscillating factors 
in the first two terms in the expressions (2) and (13) for P e2 and Piv^ — > v e ) are maximized. 
If in addition inequalities (10a) and (10b) or (14) hold, this produces the resonance maximum 
of interest in P e2 or P[y a — > u e ). For these two reasons the term "neutrino oscillation length 
resonance" was used in to denote the effect under discussion. Although this term may 
not be perfect, it underlines the physical essence of the effect. 

Conditions (9) and conditions (10a) and (10b) are shown graphically in Figs, la - If and 
in Figs. 2a - 2b, respectively, while condition (14) is shown in Figs. 3a - 3b. Figures la - If, 
2a - 2b and 3a - 3b correspond to the case of v e — i/pM mixing. We will discuss first conditions 
(10a), (10b) and (14). As was noticed in |J, a ^ small mixing angles, sin 2 29 ^ 0.10, (10a) is 
fulfilled for 

Pman ^ \^ <r ~ ^jPmanj Pman 7^ Pmani (18a) 



where 



Pc Y e c 

— yman \ I 

Pman * e 



For the standard values of p man ,c and YJ nan ' c given earlier one has r = 2.44 and (18a) 
corresponds to p T ^ an < 17.4 g/cm 3 , which is in very good agreement with our numerical 
results (Fig. 2a). Actually, as Figs. 2a - 2b show, (10a) is fulfilled for any 9 in a subregion 



n The expression for A 2v {y T — > u T ) can be obtained from eq. (1) in ref. [j]J by formally setting 
9 = tt/2 while keeping 0' m and 0'^ arbitrary, and then interchanging sm8' m (sin 6*^) and cos9' m 

(cosO- 
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of the region determined by (18a). Condition (10b) holds in two practically disconnected 
regions (Figs. 2a - 2b). In the region where (10a) is satisfied, condition (10b) reduces to the 
simple inequality 

cos(2C - A9' m + 9) < 0. (18c) 

Thus, condition (18c) determines effectively the region where both inequalities (10a) and 
(10b) are valid at sin 2 2theta ^ 0.10. Condition (18c) can be satisfied for |jj 

Pman < Pman,c ^ Pc 1 (19d) 

as well as for 

^ Pman < p m - (196) 

Only the first region, (19a), is relevant for the transitions of solar neutrinos with parameters 
corresponding to the small mixing angle MSW solution of the solar neutrino problem |I| . At 
large mixing angles the oscillation length resonance enhancement of P e 2 can take place only 
in the second region, (19b) (see further). 

In the region p r ^ an > (2r — l)p ma n, where inequality (12) holds, condition (10b) is satisfied 
at small mixing angles for 

\ 3r — 1 

" 1 j Pman < Pman g _ ~ Pman- (20) 

For r = 2.44, (20) gives 21.7 g/cm 3 < p r ^ an < 50.3 g/cm 3 - in very nice agreement with our 
numerical results (Fig. 2a). Note that, as it follows from Fig. 2a, condition (10b) is fulfilled 
in most of the region determined by (20) even at relatively large mixing angles. 
Condition (14) holds at small mixing angles for (Fig. 3a) 

Pman,c Pc> Pman T 1 Pman ■ 

The regions determined by the two inequalities (10a), (10b) and by the inequality (14) differ 
at large mixing angles. 

At small mixing angles, sin 2 29 ^ 0.10, conditions (9) can be fulfilled only for k = k' = 
JT|. For the v?, — > v e transitions in the case of v e — Vu( T ) mixing, and for the v u — > v e and 
v e — ► u h(t) transitions, taking place along the trajectories with h = 0°; 13°; 23°, one has 
AE'X' = 7r and AE"X" = vr at sin 2 29 = 0.045; 0.05; 0.08 and p r S 9.7; 9.4; 8.5 g/cm 3 , 
respectively. As can be easily seen from Figs, la - lc, Fig. 2a and Fig. 3a, all the three points 
lie rather close to, but outside the regions determined by the inequality (10b) (or (18c)) and 
by the inequality (14). Thus, for sin 2 29 < 0.040 conditions (9) are never exactly satisfied 
for the trajectories of neutrinos crossing the Earth core [Q: one has either AE'X' = tc and 
AE"X" relatively close but different from 7T, or vice versa. This can also be seen from the 
fact that at sin 2 29 ^ 0.02, one of the two conditions in (9) is equivalent to the following 
physical condition M: 



TV 



1 1 

X 7 + ~X" 



V2G F —(Y e c p c - Y e man Pman) = V2G F (N° - iVr™), (22) 



m N 



iVg and N™ an being the average core and mantle electron number densities along the chosen 
neutrino trajectory. One can easily convince oneself using eq. (8), the values of the core 
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and Earth radii R c and R 9 and the values of Y" e c , p c , Y™ an and p man that the above equality 
cannot be exactly satisfied for the trajectories of neutrinos crossing the Earth core. 

Conditions (9) are fulfilled for certain values of sin 2 29 from the interval (0.04 — 0.10), 
but in this case inequality (10b) (or (18c)) and inequality (14), which guarantee that the 
relevant neutrino transition probabilities have maxima, are not fulfilled. Nevertheless, the 
probabilities P e 2 and Piv^ — > v e ) are strongly (resonantly) enhanced even when conditions 
(9) are only approximately satisfied due to the fact that the neutrino oscillation length 
resonance is relatively wide |l[. 

As Figs, la - lc, 2a and 3a demonstrate, the point where AE'X' = tt and AE"X" = tt, 
which is closest to the regions determined by (10a), (10b) or by (14), corresponds to h = 0°. 
As a consequence, the values of P e 2 and Piy^ — > v e ) at their corresponding maxima for the 
trajectory with h = 0° and small mixing angles (sin 2 29 ^ 0.02) are bigger than max P e 2 
and max Piv^ — > v e ) for the other trajectories crossing the Earth core: max P e 2 and 
max Piv^ — * v e ) decrease with the increase of h. This is illustrated in Figs. 4a and 4b 
showing the dependence of the probability P e2 on h for p r = 10 g/cm 3 (Fig. 4a) and on h 
and p r (Fig. 4b) for sin 2 29 = 0.01. Figure 4a clearly demonstrates also that i) the neutrino 
oscillation length resonance can take place only if the neutrinos cross the Earth core, ii) this 
resonance is relatively wide in the Nadir angle variable h as well, and that iii) the resonance 
of interest is distinctly different from the MSW resonance in the mixing as the enhancement 
of P e2 due to the latter cannot and does not take place at p r = 10 g/cm 3 [0]. 

Comparing Fig. 2a (2b) and Fig. 3a (3b) one notices that the regions allowed by the 
inequalities (10b) or (18c) and (14) for values of the resonance densities from the interval 
(19a) differ somewhat even for sin 2 29 = (0.01 — 0.10). This difference is due to the presence 
of 9 in the argument of the cosine in (18c). As a result, the points at which we have 
AE'X' = tt and AE"X" = tt for the different trajectories through the Earth core are closer 
to the region where (18c) holds than to the region of validity of (14). Correspondingly, we 
can expect that for given h ^ 28°, sin 2 29 = (0.01 — 0.10), and transitions involving active 
neutrinos only, the maxima associated with the neutrino oscillation length resonance satisfy 
max P e 2 ^ max Piv^ —>■ v e ) and/or that the width of the resonance in P e2 is larger than 
that in Piv^ — > u e ). The dependence of the probability Piv^ —>■ u e ) on E/Am 2 and h for 
three values of sin 2 29, 0.005, 0.01 and 0.05, is shown graphically in Figs. 5a - 5c. 

Similar analysis can be performed for the z/2 — > u e small mixing angle (solar) neutrino 
transitions in the Earth in the case of v e — v s mixing. Equations (18a), (18c), (19a) - (19b) 
and (20) - (22) and the discussions associated with them remain valid in this case after a 
minor change: we have to replace Y e by 0.5(3V e — 1) and use for the ratio r the expression 

r = h ^ - V (23) 
— n CiY man - 1 1 K ' 

Pman\ 01 e 1 ) 

in these equations. The numerical values of the parameter r and of the upper and lower limits 
(18a) and in (20) now read, respectively: r = 2.180, p^an < 15.1 g/cm 3 and 18.3 g/cm 3 < 
PmL < 30.4 g/cm 3 . The lines of AE'X' = T\(2k + 1) and AE"X" = ix(2k' + 1) in the 
Am 2 /E — sin 2 29 plane are shown in Figs. 6a - 6c for the trajectories with h = 0°; 13°; 23°, 
while the regions where conditions (10a) and (10b) are satisfied are depicted in Fig. 7. 
As we have already indicated JT|, the oscillation length resonance conditions (9) and (10a) 
- (10b) are not even approximately fulfilled simultaneously for the — > v e transitions of 
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interest at small mixing angles. This is well illustrated by Figs. 6a - 6c and 7. Conditions 
(9) with k = k! = can both be satisfied only for h ^ 15°. At h = 16°, for example, 
we have AE'X' = tt and AE"X" = tt at sin 2 26 ~ 0.01 for Am 2 /E ~ 5.5 x 10~ 7 eV 2 /MeV 
(p r ~ 14.3 g/cm 3 ), but this point is located rather far from the region where (10a) - (10b) 
hold (Fig. 7). With the increase of h the point where one has AE'X' = tt and AE"X" = tt 
is moving further away from the region where (10a) - (10b) are fulfilled (Fig. 6c). As a 
consequence, the effect of the oscillation length resonance on the v% —>■ v e transitions at 
small sin 2 26 is considerably weaker if the latter are triggered by v e — v s mixing in vacuum. 
In spite of that, the v 2 — ► v e transition probability is noticeably enhanced in the region of 
the MSW effect in the Earth core by interference terms present in P e2 in addition to the 
MSW terms and the dominant term associated with the oscillation length resonance - 
notably by the contribution of the last term in eq. (2). Analogous results are valid for the 
v e v s transitions at small mixing. 

It is quite interesting that, in contrast to the case just considered, the oscillation length 
resonance conditions (9) with k = and k' = can both be satisfied for the Vy — > v s 
transitions at small mixing angles (sin 2 26 ^ 0.10) for the trajectories with h ^ 15°. These 
transitions are strongly enhanced for h ^ 15° in spite of the fact that the points where 
the equalities AE'X' = tt and AE"X" = tt hold lie outside the region determined by the 
inequality (14) (Figs. 8a, 8b, 9, lOa-lOb). For h = 0°, for instance, one has AE'X' = tt, 
AE"X" = tt for sin 2 26 ^ 0.05 and Am 2 / E ^ 6 x 10~ 7 eV 2 /MeV (p r = 15.8 g/cm 3 ). 
However, the dominating absolute maximum of Piv^ — > v s ) occurs in the region of the 
MSW resonance in the Earth core, i.e., at Am 2 /E ~ 4.1 x 10 -7 eV 2 /MeV (p r ~ 10.8 g/cm 3 ). 
Actually, at the point where conditions (9) hold one has Pipy — > v s ) = sin 2 (26^ — 40^J, 
but the function sin 2 (26*^ — 4:6' m ) takes a relatively small value and Pip^ — > v s ) is strongly 
suppressed. At the same time in the example we are considering max Pip^ — > u s ) = 0.94. 
The strong enhancement of Pipp, — > v s ) at small mixing angles is due to the term in the 
expression for Pip^ u s ), describibg the MSW effect in the Earth core and to an equally 
important interference term - the analog of the last term in eq. (13) for P{p i1 — > u e ). Thus, 
the mechanism of enhancement of Pipy, — > v s ) at small mixing angles is similar to the 
mechanism of enhancement of P e2 in the case of u e — v s mixing. Note that, e.g., for h = 0°, 
sin 2 26 = 0.01 and Am 2 = 5 x 10 -4 ; 10~ 3 ; 4 x 10" 3 eV 2 , the dominating absolute maximum 
in Pipy — > v s ) takes place at E = 1.2; 2.3; 9.6 GeV, which fall in the intervals of energies 
of Vy contributing either to the sub-GeV or to the multi-GeV p— like samples of Super- 
Kamiokande atmospheric neutrino events. Obviously, the discussed strong enhancement of 
Pipp, ~^ v s ) (or of Pipy — * v s ) if Am 2 < 0) can be utilized to perform a rather sensitive search 
for, and might allow to detect, Vy — > v s (or Vy — > v s ) transitions of, e.g., the atmospheric Vy 
(or Vy) at small mixing angles. For sin 2 26 = 0.01; 0.05 the dependence of the probability 
Pipp. ~> v s) on h and E/Am 2 is shown in Figs. 10a and 10b. 

At large mixing angles, sin 2 26 <; 0.10, and v e — Vyt T \ mixing, inequalities (10a) and (10b) 



12 Among the lines on the Am 2 /E — sin 2 26 plane determined by the conditions AE'X' = 7r(2fc+l) 
and AE"X" = ir(2k' + 1), those corresponding to k = and k' = are located at small mixing 
angles closer to the region where conditions (10a) and (10b) (or (14)) hold, than the lines for k ^ 
and k' ^ 0. 
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or (14) can be satisfied only in the interval (19b) (Fig. 2a). The corresponding regions of 
maxima of P e2 and Piy^ — > u e ) differ somewhat (compare Fig. 2a (2b) and Fig. 3a (3b)). 
In the regions where (10a) and (10b) hold or (14) is satisfied, conditions (9) are fulfilled 
only for k = and k' = 1, and only for h £ 20°: we have AE'X' = vr and AE"X" = 3vr 
for the trajectory with h = 23°, for instance, at sin 2 29 = 0.88 and p r = 1.23 g/cm 3 
(E/Am 2 = 3.695 x 10 6 MeV/eV 2 ). Correspondingly, the probability P{p^ — > u e ) has a 
maximum at the indicated point due to the neutrino oscillation length resonance 

Conditions (9) can be satisfied for k > and k' > 1 at large mixing angles in the region 
Pm S > as we ll) as Figs, la - lc illustrate. For the trajectory with h = 13°, for example, 
we have AE'X' = vr and AE"X" = 3tt at sin 2 29 S 0.52 and p r = 5.4 g/cm 3 ; AE'X' = 3vr 
and AE"X" = 5vr at sin 2 2# = 0.4 and p r S 15 g/cm 3 ; and AE'X' = 5vr and AE"X" = 9vr 
at sin 2 2$ = 0.42 and p. r = 23 g/cm 3 . All these points correspond, however, to minima 
in P e2 : now conditions (12) (instead of (10a)) and (10b) are fulfilled. They correspond to 
saddle points in Piv^ — > z/ e ). Some of the minima can be very deep. This is clearly seen 
in Fig. 11 showing the dependence of P e 2 on p r for sin 2 29 = 0.4. Only the minimum at 
p r ~ (6.0 — 7.0) g/cm 3 (/i = 0°; 13°) is an absolute one: we have P e 2 = at this minimum. 
Note that such a deep minimum is not present in the same region for the trajectory with 
h = 23°, in accord with the results shown in Figs, la - lc. Nevertheless, the indicated 
minimum in P e2 at h & 16° leads to the "island" of somewhat reduced values of the D-N 
asymmetry in the vicinity of the point sin 2 29 = 0.45, Am 2 = 5 x 10~ 5 eV 2 in the core sample 
of the solar neutrino events in the Super-Kamiokande detector (see Fig. 3c in ||). The wide 
deep minimum in P e2 in the case of v e — > v s transitions of the core-crossing solar neutrinos, 
taking place at sin 2 29 = (0.30 — 0.55) and p r = (6 — 10) g/cm 3 (Fig. 12), has apparently a 
similar origin Q The other minima in P e2 at p r > 10 g/cm 3 are relatively shallow. 

It should be noted that the large mixing angle (LMA) MSW u e — > transition solution 
of the solar neutrino problem is possible, as the analyzes of the solar neutrino data show, for 
Am 2 ^ 10~ 5 eV 2 and 0.60 £s sin 2 29 ^ 0.96. As we have seen, the oscillation length resonance 
maxima in P e2 can take place at large mixing angles only for p r < p ma n — 4.5 g/cm 3 . 
Correspondingly, P e2 can be enhanced due to this resonance only for neutrino energies 

E k, 30 MeV cos 29, (24) 

which can fall in the range of (5 - 14) MeV, causing substantial D-N asymmetry in the 
spectrum of 8 B neutrinos, only for 0.80 £s sin 2 29 ^ 0.95 and only for the trajectories with 
h ^ 20°. Figure 13 illustrates the dependence of the probability P e2 on p r for sin 2 29 = 0.9. 

Let us emphasize in conclusion of this Section that for most of the neutrino trajectories 
crossing the Earth core the neutrino oscillation length resonance takes place at values of the 



13 Conditions (9) (and (14) or (10a) - (10b)) are also satisfied for the v^u) — * ^etu) an d V2 — ► v e 
transitions at sin 2 28 = 1. The study of this specific case will be reported elsewhere. 

14 The existence of this minimum is reflected in the prediction of a spectacularly large in absolute 
value but negative D-N asymmetry at sin 2 26 ~ 0.50, Am 2 ~ 3.5 x 10 -6 eV 2 in the core sample of 
the the Super-Kamiokande solar neutrino induced events (see footnote 11 and Figs. 3.11, 3.12 and 
6a - 6d in §). 
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resonance density which differ from the values of the density met, according to the Earth 
models, in the Earth mantle and in the Earth core. 

3.2. Properties of the Probabilities P e2 , P{v^ — > f e ) (P(y e — > v^m)) and — > u a ) 
in the Oscillation Length Resonance Regions 

As one can easily check utilizing Figs. Id - If, 2b and 3b, practically all maxima of P e2 
and Piyfi — > v e ) associated with the neutrino oscillation length resonance, take place along 
the line AE'X' = n. For AE'X' = n(2k + 1), the probabilities P e2 and P(z/ U -> v e ) have for 
any the following form Jl|: 



P e2 = sin 2 9 + - [1 - cos AP"X"] [sin 2 (20^ - A9' m + 9)- sin 2 0] , AE'X' = n(2k + 1) 
P{^ -> i/ 6 ) = i> e - = ^ I 1 - cos AP"X"] sin 2 (2C - 40J, AE'X' = 7r(2k + 1). 



(25) 

0- 

(26) 

Expressions (25) and (26) describe P e 2 and P(u^ — > u e ) in the region of the neutrino oscilla- 
tion length resonance; they determine, in particular, the width of the resonance. Obviously, 
many of the properties of, e.g., the probability P e2 in the resonance region are determined 
by the properties of the function sin 2 (20^ — 49' + 9) in eq. (25). 

The dependence of the function sin 2 (20" t — A9' m + 9) on p r is shown graphically in Fig. 
14 for sin 2 29 = 0.005; 0.01; 0.05; 0.10; 0.50; 0.90. At small mixing angles, sin 2 29 £ 0.02, 
sin 2 (20^ - A9' m + 9) has three distinct maxima 0: at p r c es = p c where 20^ £ vr/2 and 
A9' m -9 = 5.69 £ 0.13tt < 20^, and at the points where 

2C = f, Y' (27fl) 



which is realized for 



res r*j Pman (07h\ 

Pman ~ 1±2V29 • 1 J 



At these two points 29'^ = tc. The same function has two distinct minima in the interval 

< 9 r Zn < Pc- at P Z s an = Pman (40 m = 7T, 20^ ^ tt) and at 

Pman — K f>man, (28a) 

where 



(286) 

The first minimum is rather deep, while the second one is relatively shallow 0]. At the first 
minimum sin 2 (20^ — 49' m + 9) = sin 2 0, and at the second 

min sin 2 (2^ - A9' m + 9)= sin 2 [vr - 20«((r - ft)" 1 + 2(« - l)' 1 ) + fll . (28c) 



man — 



For r = 2.44 (z/ e — v^u-) mixing) one finds k = 1.68, and correspondingly p™ n = up. 
7.57 g/cm 3 and min sin 2 (20^ — 49' m + 9) = sin 2 (7r — 13.30), which is in excellent agreement 
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with our numerical results 0. Note that the positions of both minima depend very weakly 
on 9 at small mixing angles. 

The oscillation length resonance maximum in P e2 occurs for given sin 2 29 ^ 0.05 and the 
different neutrino trajectories through the Earth core at values of p r = (8.0 — 9.7) g/cm 3 . The 
precise position of the maximum is determined by the factor 0.5(1 — cos AE"X") in eq. (25) 
and it does not coincide with the position of the maximum of the function sin 2 (26^ — 49' m + 9) 
|TJ . Actually, the maximum of P e2 falls in the interval of values of p r between the minimum 
a t Pman — K Pman — 7.57 g/cm 3 and the maximum at p r c es = p c = 11.5 g/cm 3 of the factor 
[sin 2 (20^ — A9' m + 9) — sin 2 #], where the latter rises steeply (Fig. 14). For this reason, in 
particular, the width of the neutrino oscillation length resonance in P e2 is determined at 
sin 2 29 ^ 0.05 by the expression (25) and not only by the function sin 2 (26>^ — A9' m + 9), as 
one might have expected on the basis of eq. (11). The same considerations apply to the 
position of the oscillation length resonance maximum and the width of the resonance in the 
probability P{y a — > u e ) in the indicated range of values of sin 2 29. 

The minimum of sin 2 (26^ — 49' m + 9) at p r ^ an = ^p ma n becomes shallower as sin 2 29 
increases starting, say, from the value 0.001. At sin 2 29 = 0.05 (9 = 0.11) it is still visible 
(Fig. 14), but min sin 2 (26^ — A9' m + 9) = 0.98. As sin 2 29 increases further this minimum 
disappears and the two maxima it separated, at p T ^ s an = p ma n/{^- — 2\/29) (29' m = tt/4) and 
at p r c es = p c (29!^ = 7r/2), merge to form one maximum. The latter is located practically at 
the position of the indicated minimum, i.e., at p^* n = upman, as this is clearly seen in the 
case of sin 2 29 = 0.10 in Fig. 14: at this point 29'^ — A9' m + 9 is just close to 7r/2 due to the 
relatively large value of 9. As sin 2 29 increases beyond 0.10, this maximum diminishes since 
29'^ - 49' m + 9 decreases and at sin 2 29 = 0.50, where 29'^ - 49' m + 9 = 0, it does not exist. 

The maximum of sin 2 (26^ — 49' m + 9) in the region p r ^ an ~ K Pman — 7.57 g/cm 3 for 
sin 2 29 ~ (0.05 — 0.10) is remarkably wide (Fig. 14). However, the oscillating factor 0.5(1 — 
cosA-E^X") in eqs. (25) and (26) goes through zero in the same region, thus reducing 
considerably the effect of the large width of this maximum on the probabilities P e 2, P {v^ — > 
v e ) and Pipp, — > v s ) (Figs. 5c and 10b). 

The position of the minimum of sin 2 (26 l ^ l — A9' m + 9) at p'^n — Pman, where 49' m = tt, is 
very stable with respect to the change of sin 2 29 up to sin 2 29 = 0.5. 

The maximum and the minimum of sin 2 (29'^ — 49' m + 9), which at small mixing angles 
are located at p^an — Pman/0- + 2\[29) and at p r ^ an = p ma n, respectively, "survive" at large 
mixing angles, but change their positions moving to somewhat smaller values of p r ^ n - These 
are the only extrema sin 2 (26^ — 49' m + 9) has at large mixing angles. At sin 2 29 ^ 0.4 the 
function [sin 2 (26^ — A9' m + 9) — sin 2 9] "shrinks" essentially to the region ^ pman, where 
it can take relatively large values (Fig. 14). This is also the region where the two extrema 
of [sin 2 (2#'^ — 49' m + 9) — sin 2 9] are located: the maximum occurs at p r ^ 2.0 g/cm 3 , while 
the minimum takes place at p r = (3.5 — 4.5) g/cm 3 , depending on the neutrino trajectory. 
Outside the indicated region [sin (20£ - A9' m + 9)- sin 2 9} is close to 0. 

At p r c es 3> p c and at p™ n pman we have for small mixing angles sin 2 (26*^ — 49' m + 9) = 



15 If v e — v s mixing takes place we have r = 2.18, k = 1.58 and this minimum occurs at p^an — 
7.1 g/cm 3 . In the case of — ► v s transitions, r = 2.67 and k = 1.76, so the minimum is at 
pZln = 8-0 g/cm 3 . 
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sin 2 (6 l (l — e)), where e > is a small quantity, e < 1, so that sin 2 (6 l (l — e)) = sin 2 ^. 
Nevertheless, one always has sin 2 (26^ — AB' m + 9) < sin 2 9 in the indicated regions, although 
the difference between the two sines is very small. Obviously, the probabilities P e 2 and 
P{ u n ~ * v e) are rather strongly suppressed at small mixing angles when p r c es ^> p c or p^n ^ 
Pman- The same conclusion is valid at large mixing angles, sin 2 29 ^ 0.2, for the probability 
P e 2', it is valid also in the region p r ^ s an <C pman for the probability P{up — > u e ). However, 
as can be shown (and Fig. 14 suggests), P{vp — > u e ) can take relatively large values at 
sin 2 2(9 ^ 0.2 for p J c es ^> p c . The behavior of the probability P(up — > u e ) as a function of p r 
for several different values of sin 2 29 and for h = 0°; 13°; 23° is illustrated in Figs. 15a - 15c. 

Let us note finally that the probabilities P e 2 in the case of u e —v s mixing, P e 2 jS , P{v e — > v s) 
and P(pp — > u s ) are given by somewhat more complicated expressions than (25) and (26) in 
the regions of their strong (resonance- like) enhancement at small mixing angles [|TJ: 

PIZ = \ [1 - cos AE"X"] sin 2 (2C - 0) 

+ - [1 - cosAE"X"] [1 - cos AE'X'] [sin 2 (2C - + 9)- sin 2 (2^ - t 

+ - [cos(AE'X' - AE"X") - cos{AE'X' + AE"X")\ sm(A9' m -29) sin(29 m -29'J. (29) 

The expressions for P(u e — > u s ) and -P(^ t — > v s ) can be obtained by setting formally 9 = 0, 
but keeping 9' m ^ and 6*^ 7^ in eq. (29) and using the definitions of p r ^ an c given at the 
end of Section 2. The first and the last terms in eq. (29) and the corresponding terms in 
the expressions for P{y e — > v s ) and P{v a — > v s ) give the dominant contributions. Obviously, 
they determine the properties of P e 2, s , P{y e — > z/ s ) and P(^ M — > v s ) in the region of interest. 

4. Conclusions 

The neutrino oscillation length resonance should be present in the V2 — > z^e transitions 
taking place when the solar neutrinos cross the Earth core on the way to the detector, if 
the solar neutrino problem is due to small mixing angle MSW u e — > transitions in the 
Sun. The same resonance should be operative also in the — > v e [y e — >■ f M ( T )) small mixing 
angle transitions of the atmospheric neutrinos crossing the Earth core if the atmospheric 
Up and Up indeed take part in large mixing vacuum u^ <-> u T , v n v T oscillations with 
Am 2 ~ (5 x 10~ 4 — 5 x 10 -3 ) eV 2 , as is strongly suggested by the Super-Kamiokande atmo- 
spheric neutrino data, and if all three flavour neutrinos are mixed in vacuum. The existence 
of three-flavour-neutrino mixing in vacuum is a very natural possibility in view of the present 
experimental evidences for oscillations/transitions of the flavour neutrinos. In both cases the 
oscillation length resonance produces a strong enhancement of the corresponding transitions 
probabilities, making the transitions observable even at rather small mixing angles M. Ac- 
tually, the resonance may have already manifested itself in the excess of e-like events at 
— 1 < cos9 z < —0.6 observed in the Super-Kamiokande atmospheric neutrino data [l],[n].|13[ . 



And it can be responsible for at least part of the strong zenith angle dependence present in 
the Super-Kamiokande multi-GeV and sub-GeV p— like data p] JTT| . |r3"| . 



In view of the important role the oscillation length resonance can play in the interpre- 
tation of the results of the experiments with solar and atmospheric neutrinos as well as in 
obtaining information about possible small mixing angle oscillations/transitions of the (at- 
mospheric) Up (up) and u e (z/ e ), we have performed in the present article a rather detailed 
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study of the conditions under which the resonance takes placs in the v% — *> v e , and in the 
V/j, — > v e [y e — > Vn( T )) and Vy — > v s transitions of neutrinos in the Earth. We have examined 
also some of the properties of the resonance. This was done under the simplest assumption 
of existence of two-neutrino mixing (with nonzero-mass neutrinos) in vacuum: v e — Vy( T ) or 
u e — u s , or else Vy — u s . However, in many cases of practical interest the probabilities of 
the relevant three- (or four-) neutrino mixing transitions in the Earth reduce effectively to 
two-neutrino transition probabilities for which our results are valid. Thus, our findings have 
a wider application than it may seem. 

Using our Figs, la - If, 2a - 3b, 6a - 8c and 9, in particular, one can determine practically 
all relevant values of the neutrino oscillation parameters, E, Am 2 and sin 2 26, at which the 
oscillation length resonance occurs for a given neutrino trajectory through the Earth core 
in the cases of V2 v e , Vy — > v e [y e — > Vy< T \) and Vy — > v s transitions. The sets of such 
values for the transitions of solar and atmospheric neutrinos can be easily identified. We 
have shown, for instance, that the oscillation length resonance conditions (9) with k = and 
k' = can both be satisfied at small sin 2 29 for the z/ 2 — > u e transitions in the case of v e — v s 
mixing only for h ^ 15° (Figs. 6a - 6c). In contrast, conditions (9) can be simultaneously 
fulfilled at sin 2 29 ^ 0.10 for the Vy — ► u s transitions at h ^ 15° (Figs. 8a - 8c). However, the 
points where conditions (9) hold lie outside the regions determined by the inequalities (10a) 
- (10b) and by (14). Nevertheless, as like P e 2, s 0) the probability of the Vy — ► v s transitions 
is strongly enhanced at small mixing angles in the region of the MSW resonance in the Earth 
core (Figs. 10a - 10b). This resonance- like enhancement is due to the interplay of certain 
interference terms present in the expression for Pipy — > v s ) - notably by the analog of the 
last term in eq. (13), and the term describing the MSW effect in the Earth core. For, e.g., 
h = 0°, sin 2 26> = 0.01 and Am 2 = 5 x 10" 4 ; 10" 3 ; 4 x 10~ 3 eV 2 , the dominating absolute 
maximum in -P(z^ — * v s ) takes place at E = 1.2; 2.3; 9.6 GeV, which fall in the intervals 
of energies of contributing either to the sub-GeV or to the multi-GeV fi— like samples of 
Super-Kamiokande atmospheric neutrino events. The strong (resonance-like) enhancement 
of Pipy, — > v s ) (or of Pivy — > v s ) if Am 2 < 0) found in this study can make possible the 
detection of the — > v s (or — > v s ) transitions of, e.g., the atmospheric v ^ (or v^) at small 
mixing angles. 

We have also investigated the behavior of the probabilities P e 2 and -P(^ — > f e ) (P(v e — > 
Vy( T )) in the region of the neutrino oscillation length resonance. Simple analytic expressions 
which describe P e 2 and P(z^ — * v e ) {P{ v e —* ^(t)) i n the region of the resonance and 
determine, in particular, the width of the resonance, and the probabilities P e 2, s , Pipe — ¥ u s 
and Pipy — > p s ) i n the region of their resonance-like enhancement at small mixing angles, are 
identified. We have shown that at small mixing angles i) the resonance is relatively wide not 
only in the variable E/Am 2 , but also in the Nadir angle variable h, and that ii) the values 
of P e 2 and P(i/y —>■ v e ) {P{ v e ~^ v p0r)) a t the maxima associated with the resonance decrease 
monotonically with the increase of h. Thus, at small mixing angles the enhancement of these 
probabilities due to the neutrino oscillation length resonance is maximal for the Earth center 
crossing trajectories, i.e., for h — 0°. These conclusions are valid for the enhancement of the 
probabilities P e 2,s, Pipe v s) and Pipy — > v s ) as well. 

As the results obtained in refs. |1],|2| and the present study indicate, it is not excluded that 
some of the current or future high statistics solar and/or atmospheric neutrino experiments 
will be able to observe directly the neutrino oscillation length resonance. Experiments located 
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at lower geographical latitudes than the existing ones or those under construction |19| would 
be better suited for this purpose. An atmospheric neutrino detector capable of reconstructing 
with relatively good precision, e.g.,, the energy and the direction of the momentum of the 
incident neutrino in each neutrino- induced event, would allow to study in detail the effects of 
the oscillation length resonance in the transitions/oscillations of the atmospheric neutrinos 
crossing the Earth. 
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Figure Captions 



Figures la - If. Lines in the sin 2 29 — p r (a - c) and sin 2 29 — Am 2 /E (d - f) planes on which 
conditions AE'X' = ix{2k +1), k = 0, 1, 2, ... (dotted lines), and AE"X" = n(2k' +1), k' — 
0, 1, 2, ... (black solid lines), are fulfilled for h = 0° (center crossing, (a),(d)), h — 13° (winter 
solstice for the Super-Kamiokande detector, (b),(e)), and h = 23° ((c),(f)). The parameters 
p r and Am 2 /E (the vertical axes in figures (a) - (c) and (d) - (f), respectively) are in units 
of g/cm 3 and 10~ 7 eV 2 /MeV. The results shown correspond to u e — > u^ T ) and v 2 — > f e 
transitions in the case of u e — v^r) mixing. 

Figures 2a - 2b. Regions of values of sin 2 26 l and p r (a) or Am 2 / E (b) of validity of 
conditions (10a) (A + B), (10b) (A + C), and of both (10a) and (10b) (A) for the v 2 -> v e 
transitions in the case of u e — v^t) mixing. Conditions (10b) and (12) hold simultaneously 
in the region C. 

Figures 3a - 3b. Regions of validity of condition (14) in the sin 2 29 — p r (a) and sin 2 29 — 
Am 2 / E (b) planes (A) for the v e — > v^m transitions. 

Figures 4a - 4b. The probability P e2 as a function of the Nadir angle h ((a), upper frame), 
and of p r and h (b), for sin 2 29 = 0.01 in the case of u e — * transitions of solar neutrinos. 
The lower frame in figure (a) shows the dependence on h of the derivative of P e2 with respect 
of h in the Stacey 1977 model of the Earth. 

Figures 5a - 5b. The probability P[y e — > iV T )) = Piy^ — > v e ) as a function of the Nadir 
angle h and E/ Am 2 , for sin 2 29 = 0.005 (a), 0.01 (b) and 0.05 (c). 

Figures 6a - 6c. Lines of validity of the conditions AE'X' = 7r(2fc+l), k — 0, 1, 2, ... (dotted 
lines), and AE"X" = n(2k' + 1), k! = 0, 1,2, ... (black solid lines) in the sin 2 29 - Am 2 / E 
plane in the case of u e — > v s and u 2 — > f e transitions (z/ e — z/ s mixing). The three figures 
correspond to h = 0° (a), /i = 13° (b), and h = 23° (c). 

Figure 7. The same as in figure 2b for v 2 — > f e transitions induced by u e — v s mixing in 
vacuum. 

Figures 8a - 8c. Lines of validity of the conditions AE'X' = 7r(2/c+l), k — 0, 1, 2, ... (dotted 
lines), and AE"X" = n{2k' + 1), k! = 0, 1,2, ... (black solid lines) in the sin 2 29 - Am 2 / E 
plane for — > v s transitions and h — 0° (a), h — 13° (b), and ft, = 23° (c). 

Figure 9. The same as in figure 3b for the case of P M — > z/ s transitions. 

Figures 10a - 10b. The probability P(P^ — > P s ) as a function of the Nadir angle h and 
E/Am 2 , for sin 2 2^ = 0.01 (a) and 0.05 (b). 

Figure 11. The probability P e2 as a function of p r in the case of u e — v^t) mixing with 
sin 2 29 = 0.4. 

Figure 12. The probability P e2 as a function of p r for sin 2 29 = 0.5 in the case of u e — u s 
mixing. 

Figure 13. The same as in figure 11 for sin 2 29 = 0.9. 

Figure 14. The dependence of the function sin 2 (2^ — A9' m + 9) on p r for sin 2 29 = 
0.005; 0.01; 0.05; 0.10 ; 0.50; 0.90. 

Figures 15a - 15c. The dependence of the probability Piy^ — > v e ) = P{v e — > f M ( T )) = P ef i 
on p r for h = 0° (a); 13° (6); 23° (c) and sin 2 29 = 0.005; 0.01; 0.050; 0.10; 0.50 ; 0.70. 
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